Dynamical real-space renormalization group calculations with a new clustering 

scheme on random networks 
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We have defined a new type of clustering scheme preserving the connectivity of the nodes in 
network ignored by the conventional Migdal-Kadanoff bond moving process. Our new clustering 
scheme performs much better for correlation length and dynamical critical exponents in high dimen- 
sions, where the conventional Migdal-Kadanoff bond moving scheme breaks down. In two and three 
dimensions we find the dynamical critical exponents for the kinetic Ising Model to be z = 2.13 and 
z = 2.09, respectively at pure Ising fixed point. These values are in very good agreement with re- 
cent Monte Carlo results. We investigate the phase diagram and the critical behaviour for randomly 
bond diluted lattices in d=2 and 3, in the light of this new transformation. We also provide exact 
correlation exponent and dynamical critical exponent values on hierarchical lattices with power-law 
degree distributions, both in the pure and random cases. 

PACS. 05.50.+q, 64.60.-i 
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I. INTRODUCTION 



We have generalized the dynamical real-space renor- 
malization group (RSRG) calculations for the kinetic 
Ising model l] to dilute lattices with arbitrary number 
of nearest neighbors, motivated by an interest in the re- 
laxation behaviour on networks 2J with scale free and 
exponential degree distributions. 

We first computed the dynamical critical exponent z, 
in the Migdal-Kadanoff bond moving scheme 0, 3 on 
networks with arbitrarily high, but uniform, connectiv- 
ity and found as have others [fj, Q , that z gradually con- 
verges to unity, as the spatial dimension of the system 
becomes very large (for d = 12, z — 1 = 10~ 5 ), and the 
correlation length exponent converges to v = 1. (See Fig. 
HandEJ) This is in contrast to the expected Mean Field 
values of v = 0.5 and z = 2 above the upper critical di- 
mension. The dynamic real-space renormalization group 
calculation thus yields neither a sharp crossover above 
the upper critical dimension, d c —A, nor the correct Mean 
Field behaviour in the high dimension limit. We have 
observed that a static RSRG calculation with Migdal- 
Kadanoff bond moving scheme (MK) also converges to 
the same limits, but it does so from above, whereas the 
dynamical RSRG calculation does so from below. (See 
Fig. 13 

The reason why the MK approach fails to provide a 
reasonable approximation to the critical behaviour of d 
dimensional hypercubic lattices in large d, is twofold. 
The first is because it underestimates the contribution 
from the loops in the d-dimensional Euclidean lattice 
and this effect leads to more and more inaccurate results 
for large dimensions. We have found a way to improve 
the performance of the RS approach in relation to d- 
dimensional hypercubic lattices, by defining a new type 
of RSRG cluster, which retains the inter-connectivity of 
the moved spins, and were able to obtain a convergence 
to v = 0.63 (better than 1 but still larger than 0.5) in 



10 12 



FIG. 1: The dynamical critical exponent z versus dimension d 
obtained via the conventional Migdal-Kadanoff bond moving 
scheme, z converges to 1 for large d. 



the limit of large dimensionality (d > 13). ( See Fig. |3J) 
Likewise, the dynamical critical exponent in this scheme 
converges to the value 1.6 for dimensions d > 11. (See 
Fig. The new clustering scheme for the RSRG also 
performs very well in low dimensions. The dynamical ex- 
ponent calculated within our scheme for <i = 2isz = 2.13 
and for d = 3 is z = 2.09, to be com par ed with the best 
simulation results 0, H, H EQH 13 El , which give val- 
ues between 2.11 and 2.24 for d = 2 and between 2.01 
and 2.11 for d = 3. 

The second reason for the divergence of MK results 
from those on Euclidean hypercubic lattices, is the 
nonuniform topology of the underlying hierarchical lat- 
tice, on which the MK scheme for RSRG is realised as 
an exact transformation. [T3. Il5j We turn this feature to 
an advantage in investigating the dynamical behaviour 
of the Ising model on networks with power-law degree 
distributions. 
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FIG. 2: The correlation critical exponent v versus dimension 
d via the conventional Migdal-Kadanoff bond moving scheme. 
v converges to 1 for large d. 
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FIG. 3: The correlation critical exponent v versus dimension 
d obtained by our new clustering scheme, v converges to 0.63 
for large d (compare it with Fig. |5J. 
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FIG. 4: The dynamical critical exponent z versus dimension 
d obtained by our new clustering scheme, z converges to 1.6 
for large d (compare it with Fig. 0. 



Going over to the dynamical critical behaviour for 
bond diluted lattices, we first applied the dynamic RSRG 
technique to the bond diluted kinetic Ising system in two 
and three dimensions. The approximation we used to 
compute the configuration averages is in fact exact for 
annealed randomness. The approximate transformations 
for the dilution parameter p and the coupling constant 
are well behaved as one approaches the separatrix, or 
critical line, from above, but the approximation breaks 
down in the ordered phase and for low temperatures in 
the disordered phase. In both two and three dimensions 
we were able to compute the RG flows on the disordered 
side of the separatrix and on the separatrix itself, and 
thereby determine the phase diagrams. In two dimen- 
sions we find a disorder fixed point, which is unstable, 
and which we interpret as a tricritical point, since there 
the second order phase transition line gives way to a first 
order transition line. The pure Ising fixed point is sta- 
ble, and determines the exponents along the second order 
transition line. In three dimensions, no disordered criti- 
cal fixed point is found. The critical line is depressed to 
zero temperature at a concentration p e > p* , where p* is 
the percolation fixed point. The flow along the critical 
line is to the pure Ising fixed point at p — 1 , and thus the 
critical exponents along the critical line are the same as 
the pure Ising exponents, also in three dimensions. Com- 
puting the effective critical exponents along the critical 
line, we find that z e ff varies non-monotonically as a func- 
tion of p, within the intervals [2.01, 2.25] for d = 2 and 
[2.09,2.69] for d = 3. 

Our scheme as well as the conventional equ ilibrium 
Migdal-Kadanoff RSRG (see also Ref. [H H3) fails to 
predict the crossover to a disorder critical fixed point for 
d = 3, both demonstrated by means of finite size scal- 
ing arguments |Tsj applied to large Monte Carlo simula- 
tions ln| and expected on the basis of the Harris crite- 
rion [20j. The value we find for the pure system specific 
heat exponent a via the hyper-scaling realtion 2 — a = dv 
in d = 3 is negative for the static RSRG calculation for v, 
while the dynamic calculation yields a positive a. These 
results have to be interpreted in the context of the still 
ongoing debate on the criteria for the stability of the 
pure-system critical behaviour. The rather extensive lit- 
erature on the Harris criterion 12011 has been recently re- 
viewed by Janke and Weigel |2l|. It has been shown 
by various authors [HI l22ll23| that the Harris criterion, 
which equates the crossover exponent for randomness, <f> 
to the pure system a, is simply not applicable on hi- 
erarchical lattices, and various alternative criteria, such 
as the "wandering exponent" 22] for correlations in the 
non-periodic variations in the number of bonds incident 
on lattice points (the degree of the node) have been pro- 
posed. The calculation of this exponent for our present 
RSRG scheme goes beyond the scope of this paper and 
will be considered in a separate publication. 

The paper is organized as follows. In the next section 
we setup the dynamical RSRG calculations for bond di- 
luted hypercubic lattices, and introduce a new clustering 
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FIG. 5: The "central" spin on which we concentrate is denoted 
by Co, its nearest neighbor spins by , and the interaction 
constants by where i takes values from 1 to To . 



scheme on which we will implement it. The last section 
includes our results and a discussion of the relevance of 
our results to non-uniform lattices with power law degree 
distributions. 



jm (t) = -m„(t) + (tanh (j^K^a^) 



(3) 



where = I3j[ i] and (3 = {ksT)- 1 . Here the brackets 
( ) denote both the thermal and the configuration av- 
erages. If we let r be the number of nearest neighbors 
interacting with a spin (To, then the distribution of T is 
given by 

P ( r ) = E ( T f) ^ P)> r °~^( r - ( r o - «)) ■ (4) 



Now if we expand the function appearing inside the 
brackets in Eq. @, in terms of spin products, and then 
sum over all possible configurations, we obtain 



r=i ^ ' 



l + Jt )mo(t) = 



II. THE DYNAMICAL RSRG CALCULATIONS 
FOR BOND DILUTED HYPERCUBIC LATTICES 



£m« +g(p,K,t) 



(5) 



In order to investigate the effect of the underlying lat- 
tice of arbitrarily high degree on dynamical behaviour of 
an interacting system living on this lattice, we consider 
an Ising model on the nodes of a hypercubic lattice of 
rf-dimensions, which will be subjected to bond dilution 
to yield a random network with a Poisson degree distri- 
bution. 

The hamiltonian of the system is given by 



(i) 



where Jy is the interaction between two nearest neighbor 
spins and cr; is the spin variable which can take the values 
+1 and —1. The sum is taken over all nearest neighbor 
pairs. In a ci-dimensional hypercubic lattice the maxi- 
mum number of nearest neighbor spins with a nonzero 
interaction is To = 2d. Now we will concentrate on a 
spin and its neighborhood: we will denote this spin by 
Co, the nearest neighbor spins and the interaction con- 
stants by a[ 1 ^ and , respectively, where i takes values 

from 1 to To as shown in Fig. |5] Here J-J; is subjected 
to the distribution, 



P(J®) = P 6(J?> - J) + (1 -p)6(J?>) 
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(2) 



A. Equation of motion for the magnetization 



Here ap (K) is the coefficient of the fist order term coming 
from the expansion of the function in terms of products of 
spin variables, for a particular configuration in which spin 
00 has r interacting nearest neighbors. These coefficients 
are given by 



«rW=pil E [I )(T-2n)tanh[(r-2n)Jf] , 

n=0 ^ ' 

(6) 

where n max takes values T/2 for even T or (T — l)/2 in 
the case of odd T. The other term in Eq. (JSJ, g(p, K, t) 
comes from the higher order spin products. We see that 
we can write the equation of motion for a spin variable 
as 



(7) 



where we ignore the higher order terms and define a(p, K) 



as 



a(p,K) = Y j 



r=i 



r - 1 



p'(l- P y°- l a r (K). (8) 



If we take the Laplace transform of the equation @ we 
obtain the equation of motion for the magnetization of a 
spin as 



Using Glauber dynamics [24| we may write down the 
equation of motion for the magnetization mo = (ctq) and 



(1 + s) too [s] = a(p, K) ^ m\ 



(i) 



(9) 
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FIG. 6: The cluster in our new scheme which will, in the 
coarse grained lattice, go to the bond connecting the central 
spin to its nearest neighbor in the ith direction. The inter- 
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mediate ("middle") spins are denoted by ty 
This approach preserves the inter-connectivity of the "mid- 
dle" spins and is thus able to estimate better the number of 
different paths, or loops, contributing to the spin correlations. 
(Our figure is in color in the on-line version.) 



B. New clustering scheme 



In this subsection we will introduce a new type of 
clustering scheme other than the conventional Migdal- 
Kadanoff bond moving approach. Since Migdal-Kadanoff 
bond moving drastically underestimates the number of 
loops in the system in high dimensions, we need to in- 
troduce a new scheme preserving the inter-connectivity 
of the nodes which will be eventually decimated (see Fig. 
EJ). We will implement the transformation for the scaling 
parameter A = 2. 

According to our new scheme, in the ith direction, 

(i) 

i = 1, ...,ro, between the "corner" spins ag and <j\ 
there exists a cluster containing A = A 1 "middle" spins 
denoted by Wj, with j = 1,2, ...,N, as illustrated in 
Fig. HJ1 All the "middle" spins can be connected to the 
"corner" spins and to the rest of the "middle" spins in 
the cluster. The "corner" spins will remain after decima- 
tion, and may have a maximum number To.c = FoA of 
nearest neighbors. The "middle" spins will be decimated 
through the last step of the dynamical RG calculations, 
and the maximum number of nearest neighbors they may 
have is r 0im = A + 1. We postulate that the bond oc- 
cupation probability p, and the bond strength J remain 
invariant after the application of our new type of bond 
moving scheme, and the interaction between any pair of 
spins on the cluster obeys the same distribution as in Eq. 
(2). The distribution of the number of nearest neighbors 



is given by 

n=Q 



r' 



(l-p) n p r ^- n S(T'-(T' max -n)), 

where T' max is To.c for the "corner" spins and Fo !m for 
the "middle" spins. 

The equation of motion for the expectation value of 
the jth "middle" spin in the ith direction now becomes, 



'+^)<*f>=7(p,*) 



mo + 



where ~f(p, K) comes from the configuration average 
given by 

T ,n 



7(p,*o = E 



r'=i 



1 



r - i 



p r (l- P y<>™~ r 'a r ,(K) 



If we write down the equation of motion for the expecta- 
tion value of the "corner" spins we obtain 

1 + | mo(t) = A(p, K) X 5>f ) , (13) 

' <=1 3=1 

where A(p, K) comes from the configuration average, 

a(p,k) = J2 f r r °; c " V'(i -P) ro -- I V(JO • 



r'=i 



(14) 



C. Obtaining the dynamical RSRG equations 

Now we are ready to perform the decimation. The aim 
is to rewrite the equation of motion for the "corner" spins 
in terms of their Ath neighbors. In our case, with A = 2, 
this means obtaining the equation of motion for mo in 
terms of the m^'s. For this purpose we will write down 

a) 

all the equations obtained for the "middle" spins tTj in 
the ith direction and sum all the equations. We get 

£<-f > - 

3=1 



l + --(A-l) 7 (p,A) 



Nj(p,K) [■ 



mo + m 2 



(15) 



Thus we obtain the equation of motion for the "middle" 

spins in the ith direction in terms of the "corner" spin 

(%) 

averages mo and m 2 . If we multiply the equation of mo- 
tion C3J for m by [l + ± - (A - l)j(p, K)] we obtain 



l + --(A-l) 7 (p,A) 



l + S ]mo(t) 



A(p,K) 



r„ 

E 

i=l L 



1 



d 
dt 



(N — l)7(p, K) 



N 
3=1 
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Using Eq. (|T5)l we get 



d 



(iV-l) 7 (p,if) 



dt 



m (t) 



A(p,K)J^N 7 (p,K) 



i (») 



(17) 



and 



NT A(p, if) 7 (p, K)}m (t) 



NA(p,K) 7 (p,K)J^; 



(18) 



If we ignore the second derivative term we obtain 
{[1 - (N - l) 7 (p, if) - NT A(p, if) 7 (p, if)] + 
[2-(N-l) 1 (p,K)]j t }m (t) = 



NA(p,K)j(p,K)^r>, 



(19) 



i=i 



Now let us write this equation in a familiar form, 

ri , 2-(iV-l) 7 (p,if) rf = 

1 1 - 7 (p, if) [N - 1 + NT A(p, if)] dV 

r 

m 2 



iVAfeifhfeif) 



i - 7 ( P , if ) [jv — i + #r A(p, ^)] 

Taking the Laplace transform we get 
2-(iV-l) 7 (p, if) 



(20) 



{1 



1-7(P. *Q - 1 + NT Q A(p, if)] 



s}to [s] = 



7VA(p, if ) 7 (p, if) 



1 - 7 (p, if) [iV — 1 + AToACp, if)] ^ ' 



(21) 



We see that the equation of motion (|21J) is in the same 
form as Eq. @. We identify the second term in the 
curly brackets as the renormalized Laplace variable s. 
The coefficient in front of the summation appearing on 
the right-hand side, we identify as the coefficient a(p, if), 
expressed in terms of the renormalized variables p, if. 
We thus obtain the RG equation for the time from 



2-(iV-l) 7 (p,if) 



s l-j(p,K)[N -l + NT A(p,K)\ 
and the implicit RG transformation for if from 
NA(p,K)y(p,K) 



(22) 



a(p,K) 



1 - 7 (p, K)[N-1 + NT A(p, if)] 



R(p,K) , 
(23) 

where a(p, if ) is given by Eq. iJSJ. The transformation 
for the renormalized occupation probability p is found 



by calculating the probability f(p) of an unbroken path 
from the spin <7q to through the cluster in the ith 
direction, and is thus independent from if. Thus, the 
fixed point value for the occupation probability satisfies 



P 



f(p*) 



(24) 



Note that this implies that at each stage of the distribu- 
tion, the distribution of the bond strengths is replaced 
by a distribution of the initial binary form, Eq.Q, with 
the renormalised parameters p and if. |25j This may hide 
from view certain features of the random fixed point as- 
sociated with the full distribution. 0> E3 

The fixed point of the dynamical RG transformation 
for if is found from 



a(p*,K*)=R(p*,K*) , 



(25) 



where p*, found from Eq. (|24|l . should be substituted. 
We can evaluate the correlation critical exponent v from 



dK , 
dK 



dR{p,K) ,da(p,K) 



dK 



I- 



dK 



\p*,K 



and the dynamical critical exponent z is given by 
s 



\p*,K* 



(26) 
(27) 



III. RESULTS AND DISCUSSION 

In the foregoing we have presented the generalized dy- 
namical RSRG framework for the kinetic Ising model on 
dilted <i-dimensional lattices with a random number of 
nearest neighbors, motivated by an interest in the scal- 
ing behaviour of relaxation times on random networks. 

For the case with no bond dilution, we calculated the 
correlation length critical exponent v and dynamical crit- 
ical exponent z with our new scheme of clustering in d- 
dimensional hypercubic lattices and found that v con- 
verges to 0.63 and z converges to 1.6 for large d values 
as shown in Figs, (l.'il 111 . The numerical values are given 
in Table 1. These results are much better than those ob- 
tained by the conventional Migdal-Kadanoff bond mov- 
ing scheme which is shown in Figs. respectively. 
Note that, here, just as for the conventional Migdal- 
Kadanoff bond moving scheme, these results for the crit- 
ical point and the correlation exponent differ from those 
found directly from the fixed point of a RSRG transfor- 
mation obtained by decimating the middle spins in the 
cluster shown in Fig. We denote the latter scheme as 
the "static approach" and have reported the results for 
the correlation length exponent found in this way under 

^static ■ 

Our new scheme yields dynamical critical exponent 
values z — 2.13 and z — 2.09 in two and three dimen- 
sions, as well as the percolation exponent v v . We report 
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TABLE I: Our results for dynamical critical exponent z and 
correlation critical exponent v with respect to space dimen- 
sion d at the pure Ising fixed point. 



d 


z 


Z-^d y Yl El m i c 




2 


2.13 


0.49 


1.15 


3 


2.09 


0.51 


0.82 


4 


1.99 


0.54 


0.72 


5 


1.88 


0.57 


0.68 


6 


1.77 


0.59 


0.65 


7 


1.70 


0.61 


0.64 


8 


1.65 


0.62 


0.64 


11 


1.60 


0.63 


0.63 



TABLE II: The fixed points and the critical exponents in 
d — 2 and d — 3. The first value of p shows the pure Ising 
fixed point, the second one shows the percolation fixed point 
for each dimension d. 



d 


P 


K* 


Vp 


^dynamic 


OL dynamic 


z 


2 


1 


0.27 




0.49 


1.03 


2.13 




0.5 


0.82 


1.43 


0.48 


1.04 


2.16 


3 


1 


0.12 




0.51 


0.46 


2.09 




0.16 




1.01 









our results in Table 2. We see that the agreement be- 
tween the known value of v v = 4/3 in two dimensions 
is about as good as the result in three dimensions, with 
the best Monte Carlo values being reported as v v = 0.88 
|26j . The values of the dynamical critical exponents are 
in very good agreement with recent Monte Carlo results, 
as summarized in Table 3. 

For the bond diluted case, we first computed the RG 
flows for d = 2 and d = 3 (See Figs. 01 ■ Due to 



TABLE III: Comparison of our results with the known val- 
ues coming from different approaches. Here dynamical RSRG 
calculations are written as DRSRG and Monte Carlo studies 
are denoted by MC. 



d 


Reference 


Method 


z 


2 


Present work 


DRSRG 


2.13 




Stauffer [JJ 


MC 


2.18 




Nightingale et al. [8J 


MC 


2.17 




Li et al. [9J 


MC 


2.13 




Lauritsen et al. ^0J 


MC 


2.13 ±0.02 




Adler et al. [llj 


MC 


2.21 ±0.03 




Droz et al. [6J 


DRSRG 


1.85 


3 


Present work 


DRSRG 


2.09 




Ito et al. [12] 


MC 


2.06 




Stauffer [JJ 


MC 


2.04 




Lauritsen et al. ^Vj 


MC 


2.04 ± 0.03 




Adler et al. [11] 


MC 


2.08 ± 0.03 




Ito [13J 


MC 


2.06 ± 0.02 




Droz et al. [6J 


DRSRG 


1.45 


6 


Present work 


DRSRG 


1.77 




Droz et al. [6J 


DRSRG 


1.02 




p, 



FIG. 7: The phase diagram, K~ x = ksT/J versus p, for 
d — 2. There is an unstable fixed point at p* = 0.5, (iC) _1 = 
1.22, indicated by * in the figure. The line of points extending 
to the right of the unstable fixed point is explained in the text. 
The flow from region IV is to the high temperature, p = 
fixed point. Thus the line from (p*,T*) to the percolation 
fixed point, (p*,0) is a first order phase transition line. The 
critical behaviour on the phase boundary extending from * to 
the pure Ising fixed point at T c at p = 1, is determined by the 
latter point. (Our figure is in color in the on-line version.) 



the high temperature approximation made in the deter- 
mination of the RG transformations, these flows are well 
defined on the disordered side of the separatrix, and also 
for temperatures less than, but close to, the transition 
temperatures, but not in the whole ordered region. Nev- 
ertheless, their examination is crucial to obtain the phase 
diagrams correctly. 

For d = 2, we find that the regions I and II flow, re- 
spectively, to the disordered and ordered fixed points at 
p = 1, T — * oo and p = 1, T = 0. The flow on the sepa- 
ratrix itself is to the pure Ising fixed point indicated by 
T c on the p — 1 line. Note that the line of fixed points 
of the equation a(p,K*) — R(p,K*) extending to the 
right of (p* , K*) and coming down to zero at p e is not a 
phase boundary, although it lies close to the separatrix 
for p > p* and passes through an unstable fixed point at 
(p*, K*). For p < p* , we find that in both regions III and 
IV, the flows are to the attractive disordered fixed point 
at p = 0, T — > oo. The line connecting the unstable fixed 
point at p* — 0.5, (K*)^ 1 — 1.22, indicated by * in the 
figure, to the percolation fixed point, (p*, 0) is therefore a 
first order phase transition line, separating a region with 
finite magnetization from one with zero magnetization. 
This suggests that the unstable fixed point (p*,K*) is a 
tricritical point (TCP), with a first order phase bound- 
ary connecting this point to the percolation fixed point 
p* = 0.5 at T = 0. We have checked that along the sepa- 
ratrix, from the unstable disorder fixed point to the pure 
Ising fixed point, the magnetization is zero. (We have 
also checked that the mean field-type equation for the 



7 




P P e P* 



FIG. 8: The phase diagram, K~ x = ksT/J versus p, for 
d — 3. There is no fixed point other than the pure Ising 
one for non zero temperatures. Thus the critical behaviour 
of the system is determined by the pure Ising fixed point for 
finite T. Note that the phase boundary comes down to zero 
temperature at a concentration p e which is greater than the 
percolation fixed point, p* . (Our figure is in color in the on- 
line version.) 

equilibrium order parameter too, which one may obtain 
from Eq. © by setting all the magnetizations to be the 
same, and interpreting the brackets as purely configura- 
tion averages, gives a second order phase transition in this 
interval, with the expected value of the order parameter 
critical exponent (3 = 0.5.) The value of ^dynamic and 
z at the TCP are, 0.48 and 2.16 respectively. (Similar 
unexpected features have arisen in other phase diagrams 
obtained via RSRG treatments of systems with random 
bonds. [13 ) 

For d = 3 the dynamical RG results for v gives, via 
a = 2—dv, once again a positive value for a (although the 
static result is negative, as can be readily computed from 
the values in Table 1). However, we now find that there 
is no K which the RG relation a(p*,K*) = R(p* ,K*) is 
satisfied. Examining the flow diagram in Fig. |SJ we see 
that the phase separation line comes down to T = at 
some p c > p*, precluding such a fixed point. The flow 
in regions I and II are respectively to the disordered and 
ordered fixed points, while on the separatrix it is once 
more to the pure Ising fixed point. For very low temper- 
atures, near p e , the details of the phase boundary are not 
available, due to the same difficulty as we encountered for 
d = 2. 

Since Monte Carlo simulations are plagued by 



crossover effects along critical lines, we also computed ef- 
fective critical exponents z e ff along the critical line. For 
each given p along this line, we solved for a p-dependent 
fixed point of K under the transformation in Eq. ll'MI) 
which now becomes a(f(p),K*) = R(p,K*). We then 
substitute p and K*(p) in Eq. (|2*2*|l and evaluate z g 
from (s/s) K , = A Zctf . We find that for d = 2, z e g first 
increases from 2.13 at p = 1 till 2.25 at p — 0.75 and 
then decreases to 2.01 at p e . For d = 3 the dependence 
on p is again non-monotonic, starting from 2.09 at p = 1, 
increasing to 2.69 at p — 0.4 and decreasing to 2.54 at 

Pe- 

The calculation of the dynamical critical exponent z for 
the bond (or site) diluted quenched random Ising model 
are currently the subject of numerical studies. Recent 
Monte Carlo simulations for the dynamical critical expo- 
nent z for the bond diluted quenched random Ising model 
have only yielded an effective exponent, z e ,fr , varying be- 
tween 0.59 to 0.27 along the critical line. [28| These values 
are markedly lower than the values found here. 

For d > 4, considering either the static or the dynamic 
values of the correlation length exponent in Table 1, we 
find a negative specific heat exponents from the hyper- 
scaling relation a = 2 — dv. This suggests, although 
not conclusively [2l]], that on all these lattices, the pure 
Ising fixed point will be attractive within this approach. 
Under dilution, the critical behaviour on the second or- 
der phase boundary will be determined by the pure Ising 
fixed point. In Table 1 we also display the dynamical 
critical exponent z at the pure Ising fixed point for these 
values of d. 

It should be recalled that the exponents we have re- 
ported so far are exact on hierarchical lattices [Til Il5j 
generated by iterative replacement of each line in the 
cluster shown in Fig. H3 by the cluster itself. We may 
define an effective dimension of this hierarchical lattice 
as being given via 2 dcff ~ 1 = N, where N is the num- 
ber of intermediate spins in the cluster in Fig. [ft] and 
where we have taken the scale transformation factor to 
be A = 2. For the undiluted case, we obtain a network 
with a degree distribution which is power law such that 
n(r) - T-t, where 7 = 1 + ln[2 + (N - l)/2]/lniV. We 
may now construct random hierarchical lattices by ran- 
domly diluting each cluster, with a uniform bond occupa- 
tion probability p, to get the asymptotic degree distribu- 
tion ~ cxp[— (1 — p)r]/r! for small p. From the discussion 
above, we conclude that for effective dimension d > 4, the 
pure Ising behaviour will be observed on the critical line 
for T > 0, and p > p e {d) on these hierarchical lattices. 
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